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ABSTRACT 

Superdense massive galaxies (r e ~l kpc; M~10 n M Q ) were common in the early universe (z>1.5). 
Within some hierarchical merging scenarios, a non-negligible fraction (1-10%) of these galaxies is 
expected to survive since that epoch retaining their compactness and presenting old stellar populations 
in the present universe. Using the NYU Value-Added Galaxy Catalog from the SDSS Data Release 6 
we find only a tiny fraction of galaxies (^0.03 %) with r e <1.5 kpc and M*>8xl0 10 M in the local 
Universe (z<0.2). Surprinsingly, they are relatively young (~2 Gyr) and metal-rich ([Z/H]~ 0.2). 
The consequences of these findings within the current two competing size evolution scenarios for the 
most massive galaxies ("dry" mergers vs "puffing up" due to quasar activity) are discussed. 
Subject headings: galaxies: evolution, galaxies: Formation, galaxies: structure, galaxies: photometry, 
galaxies: fundamental parameters, galaxies: peculiar 



1. INTRODUCTION 

The discovery (Daddi et al. 2005; Trujillo et al. 2006) 
that the most massive galaxies (M*>lO n M0), indepen- 
dently of their star formation rate (Perez-Gonzalez et al. 
2008), were much more compact in the past (a factor of 
~4 at z>1.5 than their equally massive local counteparts; 
Trujillo et al. 2007, Longhetti et al. 2007, Zirm et al. 
2007, Toft et al. 2007, Cimatti et al. 2008, van Dokkum 
etal. 2008, Buitrago et al. 2008; van der Wei et al. 2008) 
forces us to face the question of how these high-z galax- 
ies have evolved into the present-day massive population. 
Most of the current theoretical work in this area sug- 
gests "dry" mergers as the dominant mechanism for the 
size and stellar mass growth of these very dense galaxies 
(Khochfar & Silk 2006; Hopkins et al. 2008). As cosmic 
time evolves, the high-z compact galaxies are thought to 
evolve into the present-day cores of the brightest clus- 
ter galaxies. However, a fraction of these objects might 
survive intact from early formation having stellar popula- 
tion with old ages. In fact, in some model renditions the 
fraction of massive objects that could last without hav- 
ing any significant transformation since z>2 could reach 
1-10% (with a space density of ~10 -4 Mpc -3 ; Hopkins et 
al. 2008). Recently, however, Fan et al. (2008) have sug- 
gested an alternative scenario where the size evolution 
is related to the quasar feedback instead of merging. In 
this model, galaxies puff up after loosing huge amounts 
of cold gas due to the quasar activity. 

If the relic superdense massive galaxies exist in the 
nearby Universe (z Ss0.2) we should be able to find sev- 
eral thousands in the Sloan Digital Sky Survey Data Re- 
lease 6 (SDSS DR6) spectroscopy survey which covers 
6750 deg 2 (or a total volume of 3.73xl0 8 Mpc 3 up to 
z=0.2). With this aim we have used the NYU Value- 
Added Galaxy Catalog (Blanton et al. 2005a) to probe 
whether there is any nearby massive (M^~10 11 M Q ) 
galaxy as compact (r e ^l kpc) as those found at high- 
z, and if so, study their stellar populations to provide 
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unique clues for understanding the nature and evolution 
of these objects. In what follows, we adopt a cosmology 
of Sl m =0.3, n A =0.7 and H =70 km s" 1 Mpc" 1 . 

2. THE DATA 

Our sample was selected using the NYU Value-Added 
Galaxy Catalog (DR6) 2 . This catalog include photomet- 
ric information for a total of ~2.65xl0 6 nearby (mostly 
below z~0.3) objects. Around l.lxlO 6 of these objects 
have spectroscopic redshift determination. In addition, 
the catalog contains information about effective radii 
(Blanton et al. 2005b) and stellar masses (Blanton & 
Roweis 2007) based on a Chabrier (2003) initial mass 
function. 

With 0<z<0.2 and M*>8x 10 10 Af Q there are 152083 
targets. From these, 253 objects (^0.17% of the above 
subsample) have r e <1.5 kpc according to the catalog. 56 
of these objects are likely to be point sources (i.e. stars) 
since their sizes estimates in the r-band are smaller than 
0.05 arcsec and are not considered in what follows. This 
left a total number of 197 targets (~0.13%). As expected 
because of their compact nature, the vast majority of 
these objects are QSOs (146) with a tiny contamination 
of stars (3) according to the SDSS classification spec- 
tra, and are not considered in our analysis since their 
sizes are not representative of their host galaxies. Af- 
ter a detailed visual inspection of the remaining set of 
objects (48; i.e. only a 0.03% of our initial subsample) 
we find that 3 are very close to very bright stars, 8 are 
in close pair with other galaxies and 8 are edge-on disk 
galaxies where dust effects can be relevant, preventing 
us from using them for a further analysis. Consequently, 
we remain with a final selection of 29 galaxies (see Table 
1). The mean redshift of these galaxies is ~0.16 with 
a r.m.s. of 0.02. These objects have a mean effective 
radius of ~1.3 kpc (r.m.s. 0.15 kpc) and a mean stellar 
mass of ~9.2xl0 10 M Q (r.m.s. 1.2x 10 10 M Q ). Assum- 
ing a spherical symmetry this implies a stellar density 
of ~5xl0 9 M©kpc -3 (this is just a factor of ^2 smaller 

2 http:/ /sdss. physics. nyu.edu/vagc-dr6/vagc0/ 
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than some globular clusters) and a stellar surface mass 
density of CT5o~9x 10 9 M Q kpc~ 2 . Fig. Q] illustrates the 
position of the selected compact galaxies in the stellar 
mass-size plane, as well as some examples of the galaxies 
in our sample. 

2.1. Robustness of stellar mass and size estimates 

To address the robustness of the stellar mass and the 
sizes of the 29 compact galaxies of our sample, we have 
remeasured those quantities using different codes than 
those employed in the NYU Value-Added Galaxy Cat- 
alog. To check the sizes we used GALFIT (Peng et 
al. 2002). To calculate r e , we used a Sersic 2D model 
convolved with the PSF of the image (obtained from 
a nearby star to the source). After circularizing the 
GALFIT sizes we find a good agreement between both 
estimates, being the GALFIT measurement marginally 
smaller: Ar e /r e =-0.17±0.14. The Sersic indices are also 
very similar using both codes: An/n=-0.04±0.21, and 
with <n>~4.7. In addition, to check whether our galax- 
ies can be affected by a potential bias due to a recent 
central burst which resulted in a understimation of our 
sizes, we have compared the sizes found in the r-band 
against the sizes measured in the other SDSS filters (u,g,i 
and z). In all the bands the galaxies show very similar 
sizes. If anything, the galaxies are slightly larger in the u- 
band (i.e. against the hypothesis of a central starburst): 
<(r e , r -r e ,ii)/re,r> = -0-19±0.07, although the result is not 
statistically significant. 

The stellar masses were remeasured using the Bell et 
al. (2003) prescription based on the restframe (g-r) color 
assuming a Kroupa IMF. We obtain a very good agree- 
ment with the NYU Value-Added Galaxy Catalog esti- 
mates: AM*/M*=0.02±0.08. In addition to the stel- 
lar masses, we have compiled the central velocity dis- 
persions of these objects provided by the SDSS archive. 
We get a median central velocity dispersion of 196 km/s 
with a standard deviation of 30 km/s. We have also 
checked these values using the Penalized Pixel-Fitting 
method developed by Cappellari & Emsellem (2004) and 
we find an excellent agreement with the SDSS measure- 
ments Act/ct=0.00±0.07. 

2.2. Control Sample 

In order to make a consistent analysis of the stel- 
lar population properties of the compact galaxies we 
have created, using the same catalog, a control sam- 
ple of galaxies with similar stellar masses but with 
sizes representative of the average galaxy population. 
These galaxies were selected using the following criteria: 
0.8<M*<1.2xl0 11 M Q , 4<r e <6 kpc and within a volume 
of radius 30 Mpc in relation to the compact galaxies (to 
assure the environmental conditions are similar). This 
produces a control sample of 299 objects after rejection 
of galaxies closer to bright stars and undergoing merg- 
ers. The control sample shows a median central velocity 
dispersion of 180 km/s with a standard deviation of 34 
km/s, and a Sersic index <n>~3.8. Interestingly, the 
mean central velocity dispersion of the control sample is 
slightly smaller than the mean central velocity disper- 
sion of the compact sample. This trend towards larger 
central velocity dispersion, at a given fixed mass (lumi- 
nosity) , for those galaxies with smaller sizes is also found 
in Bernardi et al. (2006). 



2.3. SDSS spectra 

The stellar population properties of our sample were 
analyzed using the spectra available in the SDSS archive. 
The SDSS spectra cover 3800<A<9200 A and the fiber 
has 3" in diameter. This implies that for a typical com- 
pact object of our sample we cover the inner 3r e with a 
S/N of ^30. A stacked spectrum of the galaxies in our 
sample is shown in Fig. O We also show the stacked 
spectrum of the galaxies in our control sample. For a 
typical galaxy in this control sample, the fiber covers the 
inner lr e with a S/N-20. 

2.4. AGNs 

A potential source of error in the size determination of 
any galaxy is the presence of an active galactic nucleus 
(AGN) in its center which can bias our measurements 
towards smaller sizes. To address this issue we have es- 
timated which fraction of our compact galaxies are po- 
tential AGNs. Following Kauffrnann et al. (2003) we 
have performed the BPT diagram (Baldwin et al. 1981) 
log([OIII]/H^) vs. log([NII]/H Q ) to identify the AGNs in 
our sample. Emission line fluxes of [OIII], H^, [Nil] and 
H a are measured by fitting small regions of the spectrum 
with width <500A around the desired lines. The spec- 
tral region is fitted using a linear combination of a gaus- 
sian (to model the emission line) plus ~40 stellar popula- 
tion models for the stellar absorption spectra. Based on 
the MILES stellar library (Sanchez-Blazquez et al 2006; 
Cenarro et al 2007), these models are an extension of 
those in Vazdekis (1999; hereafter V99+). Only objects 
where the four emission lines are detected with S/N>3 
are considered. We use the demarcation to separate be- 
tween AGNs and starbursts proposed by Kauffrnann et 
al. (2003; their Eq. 1). On doing this we find only 2 
galaxies in the compact sample that could be considered 
as AGNs. Consequently the fraction of AGN galaxies in 
our sample is very low (~7%) and our results are not al- 
tered by them. Following the same criteria, we find that 
the number of AGNs objects in the control sample is 29 
(i.e. ~10%), similar to the fraction found in the compact 
sample. This result reinforces the idea that the presence 
of an AGN is not affecting our size estimates. 

3. STELLAR POPULATIONS 

To investigate whether the distinct structural prop- 
erties of both samples of galaxies are linked to differ- 
ences in their stellar population properties, we have 
analyzed both their global spectra (Fig. [2|) and their 
luminosity-weighted ages and metallicities on the basis 
of key absorption line-strength indices and stellar popu- 
lation models. 

Fig. [2] shows that the average compact galaxy looks 
clearly younger than the average control galaxy. This is 
well seen from its bluer continuum and from the fact that 
its absorption line is much stronger. Consequently, 
the metal lines of the compact galaxies are weaker than 
those found in the galaxies of the control sample. A 
more detailed disentangling of age and metallicity ef- 
fects is presented in Fig. [H This figure shows the in- 
dices H^q -an optimized age indicator by Cervantes & 
Vazdekis (2008)- and [MgFe] -mainly sensitive to overall 
metallicity (Gonzalez 1993)- measured for the compact 
(filled squares) and control galaxies (dots). All galax- 
ies in the range 0.132<z<0.159, with values subject 
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to be affected by strong A5577A [Oi] skyline residuals, 
have been rejected from the figure and subsequent anal- 
ysis. To avoid systematics among the galaxy indices due 
to different velocity dispersions, all the galaxy spectra 
were previously broadened to match the largest a value 
of the whole sample, which is 320 km/s. Thus, taking 
into account the SDSS instrumental broadening, all the 
indices in Fig. [3] are given at an overall a of 332 km/s. 
Also, V99+ models are overplotted at the same spectral 
resolution of the data. 

Index uncertainties have been measured for each 
galaxy from the error spectra provided in the SDSS 
database. The open square and circle in Fig. [3] indicate, 
respectively, error- weighted mean indices for the compact 
and control galaxy samples. For each of them, error bars 
in thick lines illustrate the la typical index uncertainties. 
In turn, thin error bars provide the error-weighted r.m.s. 
standard deviations. The fact that the typical index er- 
rors are smaller than their standard deviations seems to 
indicate that, within each galaxy subsample, there exist 
certain differences -not explained by errors- among the 
stellar populations of their individual galaxies. However, 
these differences are, by far, negligible as compared to 
those between compact and control galaxies. As already 
suggested by Fig. [2l it is clear from Fig. [3] that compact 
galaxies are much younger (~ 2 Gyr) than the control 
galaxy sample, which is typically old (~ 14 Gyr). There 
are also hints for compact galaxies having higher metal- 
licities ([Z/H] > +0.2) than the control sample, with 
averaged values slightly below solar. 

It is worth noting here that the above age and metal- 
licity results are difficult to be driven by the larger aper- 
tures (a factor of 3 in units of r e ; see Section 2.3) sampled 
for the compact objects, as detailed studies of age and 
metallicity gradients in normal ellipticals (e.g. Sanchez- 
Blazquez et al. 2006) report that their outer parts are, on 
average, less metal rich and slightly older than their cen- 
tres. Therefore, if aperture issues were significantly af- 
fecting the stellar population properties of the integrated 
spectra, we would expect compact objects to be older 
and more metal poor. Also, note that the total galaxy 
luminosity within 3r e is (for a Sersic surface brightness 
profile with a n=4) only a 1.6 larger than that within 
lr e (Trujillo et al. 2001), so if the outer regions were 
the clue for the strong differences in luminosity-weighted 
age between both galaxy subpopulations, compact ob- 
jects should have unprecedented young populations in 
their outskirts. 

4. DISCUSSION 

The surprising young ages of the superdense massive 
galaxies in the nearby universe cast some doubts about 
what are the mechanisms that most massive galaxies in 
the high-redshift universe follow to reach their current 
sizes. At high redshift the most massive galaxies are sup- 
posed to be the result of gas-rich mergers resulting in a 
compact remnant (e.g. Khochfar & Silk 2006; Naab et 
al. 2007). After this, dry mergers are expected to be the 
mechanism that moves these very massive galaxies to- 
wards the current stellar mass size relation. Within this 
scheme a non-negligible fraction (1-10%) of superdense 
massive galaxies is expected to survive intact since that 
epoch (Hopkins et al. 2008) and, consequently, they are 
supposed to have old stellar populations which is at odds 



with our findings. The above survival rate, however, is 
an upper limit since only major mergers (ratio 1:3) were 
considered at estimating this number density in Hopkins 
et al. (2008). In fact, the role of minor mergers and small 
accretion can be very relevant (Naab et al. 2007) at ex- 
plaining the puffing up of galaxies. Within the scheme 
of dry merging scenario, our result could highlight the 
importance of accounting for minor merging in order to 
make robust estimations in the number density of old 
superdense massive galaxies in the present Universe. 

If our compact galaxies are not relics of the early uni- 
verse how they could be form? A possibility is that they 
come from recent gas-rich disks merging. If this is the 
case we should not find many in the present-day uni- 
verse since the current massive disks population has not 
too much gas. This could explain why the number of su- 
perdense young massive galaxies is so scarce in the local 
Universe. 

Recently, Fan et al. (2008) have suggested a puffing 
up scenario where the superdense massive galaxies in the 
early universe can grow in size without suffering merging. 
This mechanism, analog to the one proposed to explain 
the growth of globular clusters (Hills 1980), argues that 
the observed strong evolution in size is related to the 
quasar feedback, which removes huge amounts of cold 
gas from the central regions quenching the star forma- 
tion. If this mechanism took place in all massive galax- 
ies in the past, no old superdense massive galaxy should 
be observed today as it seems to be the case. However, 
against this puffing up scenario could be the stellar pop- 
ulation ages of our galaxies. According to the Fan et 
al. model, after the quenching of the star formation due 
to the gas ejection, the galaxy needs some time to reach 
its new equilibrium configuration. Following their model, 
the amount of time is about 40 dynamical times. In mas- 
sive galaxies this will be around 2 Gyr. Interestingly, this 
is the typical mean luminosity weighted age of the stellar 
population of our galaxies. However, we see that they are 
still very compact. It is not clear how this can be fully 
accomodated within the Fan et al. scenario. 

A key element to explore the size evolution process is to 
the study the evolution of the central velocity dispersion, 
o\t, of the most massive galaxies at a given halo mass. 
The two competing scenarios, "dry mergers" vs "puffing 
up" predict a very different evolution of the central ve- 
locity dispersion as cosmic time evolves. In the merging 
scenario u* is basically constant with time, at most a fac- 
tor of 1.3 higher at high-z (Hopkins et al. 2008). In the 
"puffing up" scenario the original (before the expansion 
of the object) central velocity dispersion depends on their 
formation redshift as (l+z/o™) 1 / 2 . As the object in- 
flates, <7* changes with the effective radius as cr*ocr e -1 / 2 . 
Our results shows that galaxies in the control (already 
"puffed up" in this scenario) and compact sample have 
very similar a+. If both type of objects were formed at 
high redshift our results will be in contradiction with the 
"puffing up" scenario. If, on the other hand, the nearby 
superdense massive galaxies have been formed recently, 
our results agree with that scenario since young, low- 
z galaxies in the superdense phase should have velocity 
dispersions similar to old "puffed up" galaxies, because 
the effect of the increase of the radius essentially com- 
pensates the effect of the different Zf orm . 
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To test the prediction of the puffing up scenario is key 
to know whether the nearby superdense massive galax- 
ies are genuinally young objects. We need to constrain 
whether our relatively young mean luminosity weighted 
age estimate of ~2 Gyr is representative of the whole 
stellar population of the objects or an artifact due to a 
recent burst not involving a large amount of stellar mass. 
To address this issue we have probed the Star Formation 
History (SFH) of our objects by means of STARLIGHT 
(Cid Fernandes et al. 2005) exploring which SSP-model 
combination best fits the observed spectra. Our prelim- 
inar results show that SFHs from the compact galaxies 
are systematically different to those of the control sam- 
ple. For an average compact galaxy, more than 64% of 
the luminosity comes from stellar populations younger 
than 3 Gyr, in contrast to a 7% in the control sample. 
If this were the case, it will indicate that the superdense 
massive galaxies are genuinally (structurally speaking) 
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Superdense massive galaxies sample 
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Fig. 1. — Left: Stellar mass-size distribution of the NYU Value-Added Galaxy Catalog (DR6) galaxies. The position of our sample of 
compact galaxies is shown with circles. The stars are sources that according to the SDSS spectra classification are QSOs. Following Shen 
et al. (2003), over-plotted on the observed distribution are the mean and dispersion of the distribution of the Sersic half-light radius of the 
SDSS early-type (n>2.5; solid line) and late-type (n<2.5; dashed line) galaxies as a function of the stellar mass. The gray rectangular area 
shows the region used to extract the control sample galaxies. Right: Mosaic showing 6 typical galaxies from our sample. The solid lines 
indicate the equivalent to 10 kpc size at the distance of the objects. 
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Fig. 2. — Average mean spectra of the compact galaxies of our sample (grey thick line) and the control galaxies (thin line). The individual 
spectra were previously convolved to the highest dynamical velocity dispersion (320 kms -1 ) and divided by their corresponding means, so 
they contribute with the same weight. Each mean spectrum has been normalized by its mean flux within the spectral range showed in each 
panel. The inset shows the spectral regions around H/3 and Mg lines. 
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Fig. 3. — H^p vs [MgFe] diagram for 151 galaxies of the compact (filled squares, 18) and control samples (dots, 133) for which no 
A5577 A [Ol] skyline residuals affect the index measurements. Overplotted are V99+ SSP model grids with different ages (dotted lines) and 
metallicities (dashed lines) as given in the labels. The open square and circle are error-weighted mean indices for the compact and control 
samples respectively. Thick error bars show the la index uncertainties, whereas the thin error bars correspond to the dispersions of the 
samples. 



